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ABSTRACT 



A two-equation turbulence model has been extended to be applicable for compressible flows. A com- 
pressibility correction based on modelling the dilatational terms in the Reynolds stress equations has been 
included in the model. The model is used in conjunction with the SPARK code for the computation of high- 
speed mixing layers. The observed trend of decreasing growth rate with increasing convective Mach number 
in compressible mixing layers is well predicted by the model. The predictions agree well with experimental 
data and the results from a compressible Reynolds stress model. The present model appears to be well suited 
for the study of compressible free shear flows. Preliminary results obtained for reacting mixing layers have 
also been included. 


NOMENCLATURE 


b< 

body force of species i 


turbulence model constants 

C P 

specific heat at constant pressure 

c pi 

c p for i th species 

D{j 

binary diffusion coefficient 

Dt 

thermal diffusion coefficient 

E 

total internal energy 

U 

mass fraction of species i 

H 

source vector 

H 

total enthalpy 

h 

static enthalpy 

hi 

enthalpy of species i 

K 

reference enthalpy of species i 

k 

turbulent kinetic energy 

Mi 

molecular weight of species i 

M t 

Turbulent Mach number 

ns 

number of chemical species 

Pr,Pr t 

laminar and turbulent Prandtl numbers 

P 

pressure 

Q 

heat flux vector 

R° 

universal gas constant 
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Sc,Sct laminar and turbulent Schmidt numbers 
T temperature 

t time 

At time step 

U dependent variable vector 

U velocity vector 

Vi diffusion velocity vector of species i 

to,- species production rate of species i 

Xi mole fraction of species i 

x streamwise coordinate 

Xj j th coordinate 

y transverse coordinate 

a compressibility correction coefficient 

7 ratio of specific heats 

8ij Kronecker delta 

e dissipation rate 

e 0 compressible dissipation rate 

k thermal conductivity 

A second viscosity coefficient 

p laminar viscosity 

p t turbulent viscosity 

v kinematic viscosity 

p density 

crfc, <7 e turbulence model constants 

Tij stress tensor 

<&j flux vector in 3 th direction 

Subscripts 

a high speed stream 

b low speed stream 

t turbulent quantity 


INTRODUCTION 

In recent years, the development of an airbreathing hypersonic vehicle has received considerable attention. 
This is a complex task and requires innovative research in many technical areas such as aerodynamics, 
propulsion, structures and materials. Development of viable propulsion systems for such a vehicle is being 
pursued at numerous research institutions around the country. One such program is being carried out at 
the NASA Langley Research Center in which a highly integrated, hydrogen-fueled supersonic combustion 
ramjet (scramjet) engine is envisioned to be a viable propulsion system for hypervelocity vehicles [ 1 , 2 ]. Due 
to the complex nature of such a task, numerous research programs have been initiated. One of these research 
efforts is concerned with understanding the details of the complex flow field inside the engine and evaluating 
its main features for a range of flow conditions. The work reported here is concerned with this effort in 
which computational analysis of the scramjet flow field is carried out. 

The flow field in the scramjet combustor is highly complex. It is governed by the Navier-Stokes equations 
coupled with a system of equations describing the chemical reactions that occur. The flow is expected 
to be turbulent in most part of the combustor. This requires that the analysis be capable of addressing 
compressible turbulent reacting flows. The interaction between turbulence and chemical reactions is an 
important issue in the analysis. The presence of turbulence in any flow complicates accurate analysis of 
the flow due to the wide range of length and time scales. Exact solution of complex flows, such as that 
in a scramjet combustor, is impossible at present because of these turbulence scales. Hence some form of 
simplified treatment of the turbulence is mandatory in such problems and turbulence models are used in 
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this regard. Establishing turbulence models that are suitable for flows such as that of the scramjet is very 
difficult since the effects of turbulence on the flow as a whole and on the chemical reactions in particular are 
not well understood. The model should be accurate enough in predicting the physics of the problem without 
increasing the complexity of the solution procedure any more than is necessary. 

A variety of turbulence models, which have been used for many different flow configurations, exist today. 
These models range from the simplest mixing-length or zero-equation models to the most general Reynolds 
stress closures. Also, other means of analyzing turbulent flows, such as large eddy simulation, are being 
developed. A useful review of the field was conducted by Nallasamy [3]. One class of models that is widely 
used in turbulence computations today is the two-equation model. In this model, a differential equation for 
the mean turbulent kinetic energy and another for some form of the length scale of turbulence are solved 
along with the averaged forms of the Navier-Stokes equations. These two-equation models have been used 
with some degree of success for many engineering problems. They are relatively easy to implement in a 
given solution procedure and are computationally economic in comparison with the more general Reynolds 
stress models. These models do have their limitations and their applicability to any flow problem must be 
validated before using them. One of the major restrictions of most of the turbulence models in use today is 
that they have been developed for incompressible flows. Adhoc modifications to account for compressibility 
have been made to some of these models. A discussion on this subject is given by Yang et al. [4]. A fully 
compressible turbulence model of the two-equation and higher order does not exist because the averaged 
equations for compressible flows are not easily amenable to modelling using any of the known techniques. 
However, good progress has been made in this area by many researchers and it is hoped that turbulence 
models applicable to compressible flows should soon be available. In the present analysis, a two-equation 
model of turbulence has been used with a compressibility correction derived from the Reynolds stress closure 
model of Sarkar et al. [5]. The two turbulence variables are the turbulent kinetic energy and its dissipation 
rate. The governing equations are Favre averaged [6] thus including the effects of mean density variations. 
At present, the model does not include the turbulence-chemical reactions interactions and all the model 
constants of the incompressible model [7] are retained. 

The two-equation model is tested on a spatially developing, primarily supersonic, chemically reacting 
plane mixing layer. A major portion of the chemical reactions in the scarmjet combustor occur in mixing 
layers and all the complexities introduced by fluid mechanics, combustion chemistry and the interaction 
between them are retained by the reacting mixing layer. A detailed discussion on this topic is given by 
Drummond et al. [1]. Mixing layers have been studied extensively over the years. Some of the earlier work 
in this area is found in reference [8], which resulted after the conference on free shear flows held at NASA 
Langley Research Center in 1972, and in reference [9], in which the density effects on subsonic mixing layers 
have been discussed. Reference [10] gives a good picture of recent developments in the area of turbulent 
shear flows of which the mixing layer is a subset. However, as with many other compressible flows, reliable 
experimental data in the area of reacting mixing layers is limited. Availability of reliable experimental data 
for the purpose of validation is crucial in the development of turbulence models and the lack of this data is 
one of the reasons why models applicable to reacting flows are not common today. 

In the present work, the two-equation model is applied to the high-speed plane mixing layer problem over 
a wide range of flow conditions without and with chemical reaction. The nonreacting cases are computed for 
air-air systems. Hydrogen-air combustion models are used for the reacting flows. Single-step and multiple- 
step reaction systems are used here. The two gas streams that form the mixing layer are supersonic whereas 
the convective Mach number [11] of the layer varies from subsonic to supersonic values. The computer 
code SPARK, developed at the NASA Langley Research Center [12], has been used. It solves the governing 
equations using a fourth order compact scheme. The mean flow variables such as velocity, internal energy, 
the two turbulence variables and the concentrations of various species such as hydrogen, oxygen, nitrogen 
etc., are computed. Representative results are presented and pertinent flow features are discussed. The 
predictions have been compared with those of a Reynolds stress closure [13,14] and also with available 
experimental data [11,15,16,17]. Major flow characteristics such as the growth rate of the mixing layer, 
mean velocities, turbulent kinetic energy etc., are used for comparisons. 



GOVERNING EQUATIONS 


The Navier~Siok.es equations along with equations for energy and species continuity which govern flows 
with multiple species undergoing chemical reaction have been used [18]. These governing equations are 


Continuity 


Momentum 


~ + V ■ (pU) = 0 


+ V.(pUU) = Vt + pf^fibi 

dt 


where 


T = Tij 


, dui duj . 
k dxj dxi ' 




dx k ^ 


( 1 ) 

(2) 

( 3 ) 


with repeated indices indicating summation. 
Energy 


The total energy(kinetic + internal) is chosen as the dependent variable in the energy equation, given as 


djpE) 

dt 


+ V • (pUE) = V-(r-U) - V-g + pY^fibiiU+Vi) 

i = 1 

where, neglecting radiation heat transfer, 

(XjDxi w T y 


q = -/cVT + p^ifiVi + R°T^2J2 (]^)(^-^) 

»=i »=i/=i * * J 


Species Continuity 


Also, 


djpfi ) 
dt 


+ V • (ptf/<) = «* - V • (p/ify 


E = $>,/< - £ + 

»=i p 


UfcUfc 


with repeated indices indicating summation in the kinetic energy term. 

rT 


hi = h° + f c Pi d,T » = 1,2, ..., ns 
JT r 


fi 


p = pR-TY,^ 


i=l 


Mi 


( 4 ) 

( 5 ) 

( 6 ) 

( 7 ) 

(8) 
( 9 ) 


The diffusion velocities are found by solving 

XiXi 


VXj = + 

j=i p p 3=i 

Y> XiXj ,D T j D Ti VT 
A. p A . . I /. /. I T 


( 10 ) 
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Note that if there are na chemical species, then i = 1,2, ..., (ns — 1) and (ns — 1) equations must be solved 
for the species /<. The final species mass fraction f n , can then be found by conservation of mass since 

The procedures for evaluating quantities such as the specific heats, equilibrium constant for the chemical 
reactions, etc., are described in reference [1]. Also, details regarding the chemistry models, thermodynamic 
models and the diffusion models are found in reference [l]. 


Averaged Governing Equations 


Density-weighted averaging (also known as Favre-averaging [6]) is used to derive the equations which 
describe the fully turbulent flow from the above set of equations. The dependent variables, except density 


and pressure, are written in the form, 




<f> = 4, + $" 

(11) 

/V 

where the Favre-mean <f> is defined as 




4 = Ei. 

~ p 

(12) 


In this equation, the overbar indicates conventional time-averaging. Density and pressure are split in the 
conventional sense as, 


p = p + p' and p = p + p' (13) 

The following relations exist in this form of time-averaging: 



p<l>" = 0 

V = 0 (14) 


The averaged continuity and momentum equations are 

dpUi 


d J. + 

at 


dxi 


= 0 


(15) 


dpUi A dpUiUj 

+ 


dp 

dx{ 


apu'/u'j i dT± 

dxj 


+ 




at dxj 

The total energy E can be written in terms of the total enthalpy (17) as 

E 


dxi 


H - - 


Rewriting the energy equation (4) using the above and time-averaging, 


d(pH-p) + d(pH-p )Uj 


at 


dxj 


dxj 


(PU 3 - + P u"H") 


a_Qj_ 

dxj 


+ 


dUiTi 




dxj 


(16) 


(17) 


(18) 


In the above equations the body forces contribution has been omitted. Also, Qj represents the averaged 
heat flux term (equation (5)). The species continuity equation (6) is averaged similarly to get, 


dp fi dpfiUj 

at dxj 


Wl - 


dx 


:K?/" + p/.-(v i)y) 


(19) 
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In this equation, the term w,- contains correlations which have not been included in the present model. 


The equations for the two turbulence variables, turbulent kinetic energy ( k ) and the energy dissipation 
rate (e), are derived by suitably manipulating the momentum and continuity equations (equations (2) and 
(1)) and utilizing equations (11) - (14). Defining k and e as 


pu'/u'/ 
2 P 


( 20 ) 


du" du'! 
dxj dxj 

e = - 

P 

the equation for the turbulent kinetic energy can be written as 

dpk djwi. = _ m — 

n H * 3 dxi 


( 21 ) 


dt 


dxi 


du'/ 
13 dxj 




du/ 

dxi 


- g^( pu'/u'/ U 1 / +p'u/ % - u'/Tij ) 


(22) 


Substituting for Tjy and rearranging terms, 


dpk dpkUj 

dt dxi 


—mr dU < 

dZ 7 


du'/ du'/ 


P 


dxj dxj 


dxj 


(pu/u/u'J+p'u/ 6 i3 ) 


+ 


dk . 


—ft dp du'/ 

u '-'— — + p ' — - + 
’ dxi V dxi 


(23) 


d Xj 

where 0 is a set of molecular diffusion-like terms which can be neglected in high Reynolds number flows [4]. 
The first term on the right hand side of equation (23) is the energy production term, the second one is the 
dissipation term (21), the third and fourth terms represent diffusion of the turbulent energy and the next 
two terms represent the effect of compressibility. 


The equation for the dissipation rate (e) can be derived in a similar manner. It is not included here 
because it is quite long and the current model uses a simplified form of this equation. However, the modelled 
form of this equation, which has been used in the computations, will be included later in this report. 

Modelled Equations 

The k-e model achieves closure of the equations governing the turbulent flows by invoking the Boussinesq 
approximation which relates the turbulent stresses (Reynolds stresses) to the mean strain rate. Thus, the 
Reynolds stress tensor is written as, 


- pu'/u'/ 


,dUi dUj. 2 .. 

"* + aS7> “ 5 fk S,i 


(24) 


where is the turbulent/eddy viscosity defined in terms of some characteristic length and velocity scales. 
Here the length is taken to the turbulent length scale, k?/e , and the velocity scale is assumed to be k% 
leading to the following expression for p t . 

Pt = Cp p— (25) 


The correlations between the fluctuating velocity and the scalar fluctuations are modelled in a similar 
manner using a mean gradient hypothesis and a typical model is, 


pu 'P" = f t ( ^ ) 


(26) 
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where cr 4 , is a coefficient which, normally, is a constant. For <f> = /<, = Sc t , and for the static enthalpy, 

(4 ~ h), (?(}, = Pr t . 


Using the above, the averaged governing equations can be modified using models for the unknown corre- 
lations thus deriving a closed set of equations which can be solved. In this section, these modelled equations 
will be given. The mean continuity equation (15) does not require any further modelling. The momen- 
tum equation (16) has two terms (last two on the right hand side) that require modeling. The modelled 
momentum equation is, 


dpUj 

dt dxj 


dp 

dxi 


dxj 




dUj dUj 
da u dx{ 


3 ^dx k 


+ pk) %] 


(27) 


The correlation pu"H" in the thermodynamic energy equation (18) is split into its components here as 

pu"u"u" ^ 

pu'JH" = pu'jh" + ' ' 3 - + pu'y ! Ui (28) 

£1 

The modelled energy equation then is, 

d{pH-p) d{pH -p)Uj _ d ,, p Jk ,dh . fa dk 
dt dxj dxj Pr Prt dxj dxj 

Q 

+ (^7 - p % - K ' u ") u ' ( 29 ) 

where 0 * is a coefficient that appears in the turbulent kinetic energy equation. The intermediate steps that 
lead to equation (29) are straight forward and hence are not included here. The modelled species continuity 
equation is 


dp fj dpfjUj _ — ^ d ft pk . df 

dt d Xj m dXj [{ Sc Set’ dx 3 - 


(30) 


Modelling of the turbulence terms in the equations is a major area of research by itself. Details of the 
modelling of the various terms are beyond the scope of this paper. Here the models used, along with the 
relevant references, will be given. The production term is exact in its form and does not need modelling. 
Modelling of the diffusion terms in the turbulent kinetic energy equation has been a well studied area and 
the model employed here is one of the most widely used [19]. The terms identified in the previous section 
as the compressibility terms are included in the present analysis since the flows considered are compressible. 
Recently, there has been a considerable amount of activity in the area of modelling these compressibility 
effects for various turbulence models. A short account of some of these modelling efforts can be found 
in reference [4], Strahle [20] proposed a global compressibility correction for the turbulent kinetic energy 
equation. Recently, Sarkar et al. [5] proposed a model for the compressible dissipation in terms of the 
dissipation rate of turbulent kinetic energy and the local turbulent Mach number. This has been used in 
the present analysis as the model for the two compressibility terms in the turbulent kinetic energy equation. 
The compressible dissipation model is 


e c = aM?e (31) 

Here M t is the local turbulent Mach number defined as M t 2 = 2k /a 2 where a is the local speed of sound and 
the model constant a = 1.0. The modelled turbulent kinetic energy equation is [7] 


dpk dpkUj 

dt dxj 


pu'/u- 


dUi 

dxi 


pe(l + aM t 2 ) + ^rj(p + 


Pttdk - 
<7fc ' dxj 


(32) 


Modelling of the exact equation for the dissipation rate is extremely difficult, even for incompressible 
flows, due to the lack of understanding of the various complex correlation terms that are present. There have 
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been attempts at including some form of compressibility effects in the e-equation [4] but no viable model 
has emerged so far. Hence, the incompressible form of the e-equation is used in the present analysis [7], The 
modelled form of this equation is 


dpe dpeU : 
dt 


£>£i 


(33) 


where P k is the production term in the turbulent kinetic energy equation (first term on the right hand side 
of (32). This term can be written as 


Pk 


—n-uWi 

-ri'** d: 7, 

_.l ( dUi dUj 2 2.dU k 2 

^ 5 ( s7 + air' 'te 1 1 - 


(34) 


The model constants used in the analysis are (7,4=0.09, (71=1.44, (72=1.92, 0 ^= 1 . 0, cr e =1.3, Pr= 0.72, 
Prt=1.0, Sc=0.22 and Sc t =1.0. 


Solution of the Governing Equations 


Once the governing equations and required modeling are in hand, the equations are discretized and 
integrated in space and time towards steady state solutions. The governing equations are written in vector 
form as follows. 


3U 89 j 

dt dxj 


= H 


(35) 


where U is the vector of dependent variables, are flux vectors containing convective and diffusive terms 
(repeated indices indicate summation), and H are source terms containing production/dissipation terms in 
the equations. The temporally discrete form of equation (35) can then be written as 


U" +1 = U n - At 



- H"- +1 ] 


(36) 


where n is the old time level and n + 1 is the new time level. The flux terms are written at the old time 
level because the equations are advanced in real time at the smallest fluid time scale. The source terms in 
the k- and e-equations are decoupled by suitable manipulation of the ratio e/k in the present analysis. For 
example, in the fc-equation, the dissipation term is written as, 

pe = pk (j) (37) 


The term e/k is treated as a known quantity taking its value from the previous time step in the solutions. A 
similar treatment of the source term is done for the e-equation also. These nonlinear turbulence source terms 
are treated in a pointwise-implicit manner while solving the turbulence equations. The k- and e-equations 
are written in the form, 

AH 

(!- (U" +1 - U») = - At (38) 

The jacobian can take different forms depending upon how the source terms are written. In the present 
computations, this jacobian has the following terms: 


d Jh 

dpk 


= —[-D+ - (l + aM t 2 )] 
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0 


where 


dH k 

dpe 

dH e 

dpk 

8H t 

dpe 


0 


[-Ci D + C 2 -\ 



(39) 


(40) 


As an option, the chemical source term is also written implicitly [21] to alleviate the problem of stiffness 
in the governing equations when chemistry time scales become small as compared to fluid time scales early 
in a calculation or when the system approaches chemical equilibrium. The governing equations are written 
in two-dimensional cartesian form for the solutions. It must be pointed out here that even though turbulence 
is three-dimensional in nature, the contribution of the fluctuating velocity in the third dimension has not 
been incorporated in the calculations. The discretized equations are solved by means of the elliptic solver 
SPARK [12]. A fourth-order compact scheme is used in the solution procedure. Details of the code and the 
solution procedure can be found in [12]. 

RESULTS AND DISCUSSION 


A two-dimensional, high-speed mixing layer is considered in this study. A schematic of this flow problem 
is given in figure 1. The two streams coming off the splitter plate are supersonic. However, the convective 
Mach number [11] of the mixing layer ranges from subsonic to supersonic values. The two streams are 
air for the nonreacting cases whereas the higher-speed stream is comprised of a hydrogen-nitrogen mixture 
(10% H? and 90% N?) for the reacting mixing layer cases. The inlet mean velocity is assumed to have a 
hyperbolic tangent profile, thus imitating the flow that exists downstream of the splitter plate edge. Ease 
of computations prompted this choice of initial profile for the mean velocity. The initial distributions of the 
turbulent kinetic energy and its dissipation rate are chosen to be compatible with what is observed in such 
a flow problem. A constant turbulence intensity level is used in the free stream for arriving at the initial 
distribution of turbulent kinetic energy and the dissipation rate distribution is deduced from the A:-field 
using a length scale parameter. These initial distributions are shown in figure 2. The temperature and 
pressure are assumed to be initially uniform for both the streams. The free stream turbulence intensity 
proved to be a sensitive parameter, especially for the supersonic convective Mach number cases. A constant 
value of the free stream turbulence intensity (with respect to the mean velocity) has been maintained in all 
the calculations. Hydrogen-oxygen reaction systems are considered in the reacting flow calculations since 
this solution procedure is to be applied to scramjet combustor flow studies. One-step (one reaction, four 
species) and multiple-step (eighteen reactions and nine species) reactions have been considered. Different 
sets of computations were done based on the convective Mach number of the mixing layer. Here, for a given 
free stream temperature and pressure, the slower speed stream had the same velocity in all the calculations 
changing only the speed of the faster stream in order to arrive at the required convective Mach number. 

In this section, representative results of the calculations axe presented. Comparison with available exper- 
imental data is done in a limited sense. Also, the solutions are compared with the Reynolds stress closure 
predictions carried out at NASA Langley [13,14]. Typical predictions of the mean velocity, temperature and 
turbulent kinetic energy for a convective Mach number of 1.5 are given in figure 3. The free stream Mach 
numbers of the two streams are 3.2 and 6.2. These figures show the development of the mixing layer in the 
axial direction and the corresponding changes in the flow variables. A slight shifting of the layer towards 
the lower speed stream is seen here. The mean temperature increases, from its free stream value, inside the 
mixing layer. A similar variation is observed in the turbulent kinetic energy distribution. Figure 4 shows the 
width of the mixing layer (vorticicity thickness) as a function of the streamwise distance. Two values of the 
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convective Mach number have been considered. The growth of the layer is linear after the inital development 
of the flow and in thin linear region, similarity in the mean flow characteristics can be expected. This is seen 
in figure 5 where the mean velocity and turbulent kinetic energy profiles are plotted in similarity variable 
coordinates. Here y c represents the location where the two streams have a common boundary initially. 

The major aim of the present work is to extend the widely used incompressible two-equation turbulence 
model for compressible flows. As mentioned earlier, a compressible correction has been incorporated in 
the turbulence model used here. The predicted growth rates of the mixing layer with and without the 
compressibility correction model are plotted in figure 6. For comparison, the results using a global correction 
model [20] are also included. Here Cg is defined as, 

dS u a + u b 

6 dx U a -U b 

and Cg 0 is its value for incompressible flow (assumed to be at a convective Mach number of 0.1 here). The 
effect of compressibility on mixing layer is to reduce the growth rate of the layer with increasing convective 
Mach number. This trend is seen in the figure. For the cases without the compressibility correction and 
those with the global correction the growth rate seems to reach a constant value after a convective Mach 
number of about 0.5. However, the effects of compressibility are expected to be more pronounced beyond this 
point leading one to believe that the calculations without the above correction do not address the important 
problem of compressibility well. As seen in the figure, the effect of compressibility is well predicted by 
including the compressibility correction. This has been verified by comparing the predicted growth rate with 
available experimental data, as shown in figure 7. The growth rates predicted by means of the Reynolds 
stress closure[13] are also shown in this figure. As with any turbulent flow case, the experimental data here 
show a wide scatter [11,15,16,17]. This raises the question of whether the convective Mach number is the 
sole basis for comparison between the wide range of mixing layer data available in the literature. Setting 
aside this question for the moment, the figure shows that the trend of reduced growth rate with increasing 
convective Mach number is predicted very well by the k-e model. The comparison between these predictions 
and the Langley experimental data is excellent in the supersonic range of the convective Mach number. 

One of the important aspects of a two-equation model is that it is well suited for application in engineering 
flow problems from the point of view of ease of adaptability, computational economy etc. However, the model 
has been found to be lacking in certain flow cases where the Reynolds stress closure may be appropriate. 
So, one of the main steps in the present analysis has been to compare the k-e predictions with the Reynolds 
stress closure predictions in order to evaluate its applicability to the flow problem considered. This is done 
in figure 8. Figure 8a compares the predicted mean axial velocity profiles at representative locations and 
figures 8b and 8c compare the turbulent kinetic energy and dissipation rate profiles, respectively. The 
comparisons show that the k-e predictions agree very well with those of the Reynolds stress model [14] and 
hence for high-speed mixing layer flows applications the two-equation model seems to be suitable. However, 
it must be reiterated that this conclusion does not necessarily carry over to all compressible flow problems. 

Coming back to the question of whether the convective Mach number should be the sole basis for com- 
paring similar mixing layers, calculations were done with different free stream temperatures and velocities. 
Some of the results are shown in figure 9. In this figure, the free stream pressure is the same (1 atm) for 
cases 1 to 4 and it is increased by a factor of 2.1 for case 5. The results indicate some dependency of the 
growth rate on the free stream temperature. The figure also shows the equivalence between two sets of data 
(cases 1&5 and cases 2&4) thus provoking the above question. The quantity Rg, given by 

_ pU b dS U a + U b 
6 JT dx U a -U b 

is found to be nearly identical between cases 1 and 5 and between cases 2 and 4 for a given convective Mach 
number. The streams that form the mixing layer are both air in these cases. This leads to the speculation 
that in order to have equivalence between two mixing layers, the parameter Rg must also be considered 
which is a Reynolds number like parameter for mixing layer. 
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As mentioned earlier, the computations were sensitive to the initial distributions of the turbulence field. 
However, the free stream turbulence level does not affect the predicted characteristics of the mixing layer 
significantly. The compressibility correction term is dependent upon the local turbulent Mach number as 
defined in equation(31). The magnitude of this term reaches a maximum value approximately equal to one- 
third the dissipation rate (e) inside the mixing layer. This indicates that the two-equation model without 
any correction for the effects of compressibility will be grossly in error in predicting compressible flows. 

Figures 10 - 12 show the results obtained for a reacting turbulent mixing layer using the model described 
in this paper. Hydrogen-oxygen reaction is modelled by means of a 18-reaction steps, 9-species system 
in these calculations. The free stream velocities, pressure and temperature are identical between the two 
cases. The free stream temperature is 2000 K and the pressure is 1 atm. Comparison between reacting 
and nonreacting cases are shown in figures 10 and 11. Figure 10 shows the width of the mixing layer for 
a supersonic convective Mach number. Figures 11 are the mean temperature, axial velocity and turbulent 
kinetic energy distributions. Examination of these figures indicates that the temperature distribution is 
altered, as expected, due to the heat release during the chemical reactions. The temperature reaches a peak 
value of about 2510 K inside the mixing layer. However, the mean dynamic and turbulence fields do not 
change. The width of the shear layer remains almost constant. Figure 12 shows the distributions of the 
primary species in the flow. Mass concentrations of hydrogen, oxygen and water vapor are shown. The 
extent of the reaction zone can be seen clearly here. 

CONCLUDING REMARKS 

A two-equation turbulence model (A; — e) has been modified to be suitable for addressing compressible 
flows. A compressibility correction model based on modelling the dilatational terms in the Reynolds stress 
equations has been used. A two-dimensional high-speed mixing layer is studied using the model. Comparisons 
of the predictions with available experimental data and the predictions of a compressible Reynolds stress 
closure indicate that the model is well suited for the study of such flows. The decrease in the growth rate of 
the mixing layer with increasing convective Mach number is well predicted by the model. A parameter which 
may be useful to establish the equivalence of mixing layers has been identified. Representative solutions of 
reacting high-speed mixing layers also have been given. This two-equation model is being developed for 
application to reacting, compressible flows. 
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